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We calculate the decay of r~ — ► n~"/v T in the framework of Resonance Chiral Theory (R\T). By demanding the 
high energy constraints from QCD on the related form factors, we could predict the various physical observables 
of r~ — > n^jv T without any free parameter. Our results show that for a realistic cut on the photon energy 
(around 100 MeV) this mode gives a branching ratio of roughly 0.1% that should have already been detected at 
the heavy-flavour factories. Another interesting subject we have studied based on our calculation of the decay 
r~ — ¥ ■K~yv T is the experimental background estimation of the lepton flavour violation process r~ —¥ (J,~j- We 
point out that although the description of radiation that PHOTOS provides -which has been used by BaBar 
and Belle collaborations to estimate this source of background- is in excellent agreement with the theoretical 
expectations in the low energy region in r~ — ¥ n~jv T decay, this is not the case in the high energy region, 
precisely where it is easier that this decay mimics the process r~ — ¥ M~7- 



1. Hadronic decays of the r lepton 

The decays of the r lepton into hadrons are a 
very interesting subject in Particle Physics that 
has been a very active area of research during the 
last 25 years [I]. The r is privileged because its 
mass allows it to decay hadronically while being 
a lepton, which means a clean place to scrutinize 
hadronization effects at low-energies, where QCD 
becomes nonperturbative. One can exploit this 
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advantage to a great extent through its inclusive 
decays [5] • However we will focus in what follows 
in one of its exclusive decays. 
The decay amplitude for the one meson ra- 
diative decays of the r includes an internal 
bremsstrahlung (IB) component, that is given by 
QED, and thus can be calculated unambiguously 
to any desired order in perturbation theory. In 
addition, one has the structure dependent (SD) 
part, dominated by the effects of the strong inter- 
action and the interference (INT) between both. 
Lorentz symmetry determines that there are two 
independent structures, the so-called vector (V) 
and axial- vector (A) form factors that encode our 
lack of knowledge of the precise mechanism re- 
sponsible for hadronization. 

We consider the r~(p T ) — ¥ v t (c[)it~ (p)j{k) pro- 
cess @. The kinematics of this decay is equivalent 
to that of the radiative pion decay [4]. We will 
use t := (p T — q) 2 = (k+p) 2 . In complete analogy 
to that case [5] , the matrix element for the decay 

2 The decay t~ (p T ) — ¥ u T (q)K~ (p)-f(k) , which is related 
to it via chiral symmetry, is discussed in Ref. [3] , in which 
this proceeding is based. 
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of t — > 7r ju r can be written as the sum of four 
contributions: 

M [r~ -4 f T 7r~7] = -A^/b^+A^/b^+A^^+A^a , (1) 
with 3 

iMiB T = G F V ud eF^p^(k)L" v , 

iM-IB-K = G F V ud eF^(k) \ LM> 

iMv = iG F V ud eF v (t)e^ pa e 1 ' (k)k p p a V 1 , 

iM A = G F V ud eF A {t)e v {k) [(t - mj) 9p - 2p v k ll ] 

where e is the electric charge of the positron and 
e„ is the polarization vector of the photon. Fy (t) 
and -Fa(£) are the so called 5Z) form factors. Fi- 
nally L M and L^ v are lepton currents dehned by 

= ^ T (9)7 M (l-75K(p T ), 

= ^ {qW - t.) ( !:^:% ■ 

The notation introduced for the amplitudes de- 
scribes the four kinds of contributions: M.ib t 
is the bremsstrahlung off the tau, (Figure [Ha)); 
■M.ib„ is the sum of the bremsstrahlung off the 
7r (Figure [ljb)), and the seagull diagram (Figure 
0Jc)); Aiy is the SD vector contribution (Figure 
QJd)) and Ma the SD axial- vector contribution 
(Figure (He)). Our ignorance of the exact mech- 
anism of hadronization is parametrized in terms 
of the two form factors Fj\_(t) and Fy(t). In fact, 
these form factors are the same functions of the 
momentum transfer t as those in the radiative 
pion decay, the only difference being that t now 
varies from up to M% rather than just up to 
to \. The dimensionless variables 

_ 2p T ■ k _ 2p T -p 

allow to measure the photon and pion energies in 
units of i R t ne tau rest frame. Their kinemat- 
ical boundaries are given by 

0<x< 1-rl, l- x + ^<y<l + rl,(3) 
where 

^ :=(^) 2 - 0.006 «1, (4) 

3 Definitions and conventions differ with respect to Ref. [6] 
(DF), in which we have detected [3] several typoes. 




Figure 1. Feynman diagrams for the different kinds 
of contributions to the radiative decays of the tau in- 
cluding a pion, as explained in the main text. The 
dot indicates the hadronization of the QCD currents. 
The solid square represents the SD contribution me- 
diated by the vector current and the solid triangle the 
SD contribution via the axial-vector current. 



so that the photon spectrum will extend up to ~ 
0.994-^. In the decay t — > /17 the photon energy 

would be hxed to ^ (l - - 0.996^ so 

one could expect some contamination of the decay 
t~ — > 7r~7^ r to r — > fij, with the tt misidentified 
as a /x. The invariant mass of the 7T-7 system 
will vary in the range t\M^ < t < when the 
other independent kinematical variable is taken 
as x = x{z), with z — t/M% ■ 

2. Theoretical setting 

Tau lepton decays probe QCD in its non- 
perturbative regime where standard expansions 
in powers of the coupling constant are no longer 
applicable. However, the chiral symmetry of the 
massless theory allows to build an effective field 
theory dual to QCD in the light-meson sector, 
Chiral Perturbation Theory (xPT) [7 . This will 
not be enough, since t can reach which is 
far beyond the region of applicability of xPT. 
Still, the low-energy region of semileptonic tau 
decays is well described by the xPT results [8]. 
Then, one should envisage a way of enlarging the 
range of applicability of ^PT while respecting its 
low-energy behaviour. The limit of a large num- 
ber of colors (Nc) in QCD [9] is a useful tool 
in the development of Resonance Chiral Theory 
(R\T) [10], a Lagrangian formulation including 
the resonances as dynamical degrees of freedom 
that preserves assymptotic QCD properties 



|12) . Remarkably, RxT was capable of predict- 
ing the low-energy constants (LECs) of xPT in 
terms of resonance masses and couplings. The 
phenomenological application of the theory (in 
the Nc — > oo limit) to study two- [13] and three- 
meson decays of the r [14] has been successful and 
several Green- functions [15], [16] and associated 
form factors have been studied within it over the 
years. 

The leading action of the RxT Lagrangian in- 
cludes the 0(p 2 ) xPT m the even-intrinsic par- 
ity sector and the leading 0(p A ) xPT (given by 
the Wess-Zumino-Witten term 17 ) in the odd- 
intrinsic parity sector. Higher-order pieces in 
the chiral expansion are assumed to be gener- 
ated by the integration of the resonances (this 
was checked to occur at 0(p 4 ) for the Li cou- 
plings of RxT [10]). 

Next, one adds all pieces including resonances 
(R) and chiral tensors with low enough chiral 
order (x^) to not violate high-energy condi- 
tions or to force fine-tuned cancellations among 
the relevant couplings to fulfil the short dis- 
tance constraints. In the odd-intrinsic-parity sec- 
tor, that contributes to the vector form factor, 
this amounts to include all terms of Rx^ and 
RRX {2) - For the even-intrinsic-parity operators 
that are contributing to the axial- vector form fac- 
tors, these are the terms of Rx^ 2 - Since previous 
analysis showed the relevance of the RRx^ with 
negative intrinsic parity we will consider them 
here, as well. All mentioned pieces of the La- 
grangian can be found in Ref. [3], as well as the 
expressions for the SD form factors in r —> tt"/u t 
decays. In Figs. [2] and [3] the Feynman diagrams 
contributing to these form factors are displayed. 
The thick dots represent strong vertices. 




Figure 2. Axial- vector current contributions to 
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Figure 3. Vector current contributions to r 
7r~ ' ~yv T . 



We provide the resonances with an adequate 
energy-dependent width obtained consistently 
within RxT |18j and consider only one resonance 
multiplet per set of quantum numbers (the con- 
tribution of spin-zero resonances is suppressed by 
conservation laws and the fact that there is only 
one meson in the final state). Since this process 
has not been measured yet, one should try the 
simplest possible description (in order to be pre- 
dictive) that can eventually be completed if the 
data require so. 

Next we require a Brodsky-Lepage |11] behaviour 
to the vector form factor and demand that the 
axial-vector form factor satisfies a once sub- 
tracted dispersion relation. This results in con- 
straints among the Lagrangian couplings. Note- 
worthily, the relations found in the r — » ir/K"fv T 
decays are all consistent with those obtained in 
the r decays into two and three mesons. There is 
only one relation that differs with respect to the 
study of the Vector- Vector-Pseudoscalar Green's 
function |15] for a coupling whose impact is very 
mild and the discrepancy is less than 10%. 

3. Phenomenology in n~~fi' T decays 

First we have assessed the importance of the 
model independent contributions, meaning the 
IB described by QED and the Wess-Zumino- 
Witten contribution to the vector form factor, 
that is determined from QCD. We have thus 
switched off the remaining contributions to the 
vector form factor and the whole axial-vector 
form factor in this first step. We find that for a 
cut on the photon energy of 100 MeV, it amounts 
to 0.9% of the non-radiative decay, namely a 
branching fraction of 0.1%. Next, we include also 
the model dependent contributions, obtained as 
described in Sect. [2J 
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In Fig. @] we show the differential decay width of 
the process r~ — > tt~"/i> t including all contribu- 
tions as a function of x, i.e. the photon energy 
in the tau rest frame and in Fig. [3] we display 
the SD contributions, that are enhanced near the 
endpoint region, where the process t~ — >• 7r~7^ T 
can contaminate the decay t~ — > fJ>~~/. 




X 



Figure 4. Differential decay width of the process 
r~ — > tv~"/v t including all contributions as a function 
of x. 



4. Comparison with PHOTOS 

In the previous section, we have shown that for 
a realistic cut on the photon energy (around 
100 MeV) this mode gives a branching ratio of 
roughly 0.1% that should have already been de- 
tected at the heavy-flavour factories. Notwith- 
standing, this decay mode has not been measured 
yet. This does not mean that the detection of 
soft photons at B-factories is not as good as esti- 
mated, but that the splitting of the radiative and 
non-radiative pion decay of the tau was not con- 
sidered a priority. 

According to the last report by the Belle [TO] Col- 
laboration 0, the main source of background in 
t~ — > /i~7 searches is the process t~ — > 7r~7^ T , 

4 Although the BaBar report |20l is not that detailed, we 
assume similar figures to hold in their case. 
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Figure 5. Differential decay width of the process 
r~ —5- 7r _ 7^ r including only the structure dependent 
contributions as a function of x. 



where the pion is miss-identified as a muon. This 
decay represents up to 80% of all background and 
the pion fake probability to be detected as a muon 
is 0.8%. However, the most important contribu- 
tion to this background would come from radia- 
tion off the e + e~ pair and not off the 7r or r. The 
background related to the decay t~ — > 7r~7J/ T is 
estimated using PHOTOS [21 that only incorpo- 
rates the model independent contributions from 
QED to r~ — > 7r _ 7^ T , i.e. the IB parts in Eq.([T]). 
For low photon energies, this decay is dominated 
by the IB parts and is detected as t~ — > 7r~j/ T , 
since the photon cannot be resolved. However, 
we have seen in Fig. 2] that for large photon en- 
ergies the dominant contribution is given by the 
SD terms. Indeed, it is precisely near the end- 
point of the photon decay spectrum where it is 
easier that the miss-ID happens, since in the two 
body-decay t~ — > [i~-f, E 1 ~ ~ M T /2, and 
the correction is 1 =F 0.004, while in t~ — > n~^v T 
the photon energy spectrum extends essentially 
to the same value and the correction here is given 
by 1 - r\, see Eq. (gl). 

We compare our prediction for the photon spec- 
trum in r~ — > 7r~7^ T with what is obtained with 
PHOTOS [21] and MC-TESTER [22] runned 
with fixed first order only (similar results are ob- 
tained with exponentiation on), as it is displayed 



in Fig. [SI 

In order to assess the impact of the SD effects in 
that background estimation, one needs to com- 
pare PHOTOS and our prediction for the case 
where the 7r and 7 energies reconstruct the r mass 
up to 9 MeV of difference. We do this in Fig. 
[7j where this missing mass is distributed evenly 
between the it and the 7 and only the relevant 
region near the endpoint is displayed in order to 
better appreciate the differences. The PHOTOS 
simulation was obtained with lOOMevents gener- 
ated for t~ — > 7r~ (j)v T where 2.9% corresponds 
to the radiative decay. We see that in the last six 
bins there is only one event. We can estimate - 
rather conservatively- the difference between our 
prediction and PHOTOS by taking the integral 
over these last six bins, where we will estimate 
the PHOTOS contribution by our prediction in- 
cluding only IB. This gives a ratio of ~ 5, as the 
underestimation of background due to the SD ef- 
fects we have studied. 

Nevertheless, we stress that the total decay 
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Figure 7. Comparison of our prediction with PHO- 
TOS for the differential decay width of the process 
t~ — > n~jv T including all contributions (in our case) 
as a function of x. In both cases is integrated over 
the range 1 - r I - 0.05 < y < 1 - r\ + 0.05. 
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. {E 1 = lOOMeV) = 1.1. 



Our estimations (factor of 5 more background 
due to the considered decay that obtained with 
PHOTOS) would imply that this source of back- 
ground could affect the present upper bounds 
in Refs. QU, [20]:4.5 X 1CT 8 at 90% CL, mak- 
ing them even stronger. The new version 23 
of hadronic currents in TAUOLA [24] and the 
inclusion on PHOTOS of our current in the line 
of Ref. [25] will be the ideal tools to estimate 
reliably this background. 



5. Conclusions and Outlook 



Figure 6. Comparison of our prediction with PHO- 
TOS for the differential decay width of the process 
r~ — > n~"/v T including all contributions (in our case) 
as a function of x. In both cases E w is integrated over 
all its range. 



width is mildly affected by the SD parts so that 
the global description provided by PHOTOS is 
good. For example, for the total decay width one 



We have analyzed the decay r~ — > n~~fv T . It 
turns out essential to measure the photon energy 
spectrum near the endpoint in order to assess the 
real background affecting the lepton flavour vio- 
lating decay r~ — > ^~7, and the upper bound for 
its branching ratio. Moreover, it would be an in- 
teresting tool for lepton universality tests through 
the ratio T(t~ — > n~jis T )/r(Tr~ fj,~ 7i^) . To 
our view this decay channel should no longer 
be regarded as attached to the t~ — > tt~v t 
measurement for energetic enough photons which 
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would allow experimental resolution. Therefore, 
it would become a golden mode to be searched 
for by the (super)B and tau-charm factories in 
the near future. 
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